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摘 要 : 视网膜 血管 的 形态 和 结构 在 眼科 疾病 诊断 中 起 着 非常 重要 的 作用 ， 同 时 对 于 心血 管 疾病 和 糖尿 病 的 诊断 也 有 
和 要 人 作用， 视网膜 血 管 自动 分 割 是 基础 。 针 对 视网膜 图 像 采 集 过 程 中 由 于 疾病 引起 的 图 像 光照 反射 过 强 问 题 ， 提 出 了 
一 种 修正 的 形态 学 与 Otsu 相 结合 的 无 监督 视网膜 血管 分 割 算法 。 首先 运用 形态 学 中 的 高 低 幅 变换 增强 血管 与 背景 的 对 
比 度 ， 然 后 提出 了 一 种 修正 方法 ， 消 除 部 分 视网膜 疾病 引起 的 光照 问题 ， 最 后 使 用 Otsu 阔 值 方法 分 割 血管 。 算 法 在 
DRIVE 和 STARE 视网膜 图 像 数 据 库 中 进行 了 测试 ， 实 验 结 果 表 明 ，DRIVE 数据 库 的 分 割 精度 为 0.9382，STARE 数据 
库 的 分 割 精度 为 0.9460， 算 法 的 执行 时 间 为 1.6s。 算 法 能 够 精确 地 分 割 出 视网膜 血管 ， 与 传统 的 无 监督 视网膜 血管 分 
划算 法 相 比 ， 算 法 的 分 割 精度 高 、 抗 干扰 能 力 强 。 
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Abstract: The appearance and structure of retinal vessels play an important role in the diagnoses in ophthalmic diseases, it also 
plays an important role in the diagnosis of cardiovascular disease and diabetes, which requires an algorithm to automatically 
extract the retinal vessels. To enhance the image with imbalanced local illumination caused by retinal diseases in the process of 
fundus image acquisition, this paper proposed a new unsupervised retinal vessel segmentation method with morphological and 
Otsu. First, used the combine of the top-hat transformation and the bottom-hat transformation to enhance the contrast between 
the blood vessels and its background in a retinal image. Next, presented a novel revised method to remove the problem, which 
is caused by the imbalanced local illumination in the enhanced retinal image. Finally, applied the threshold calculated by Otsu 
method to extract the retinal vessels. We evaluated the algorithm with two publicly retinal image databases DRIVE and STARE. 
The experiment results indicate that the segmentation accuracy in DRIVE database achieves 0.9382, the segmentation accuracy 
in the STARE database achieves 0.9460. The run time of our new method is 1.6 seconds. The new algorithm can accurately 
extract the retinal vessels. Compared with the traditional retinal vessel segmentation algorithm, its segmentation accuracy and 
anti-perturbation ability Improve. 
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疾病 包括 高 血压 、 心 血管 病 、 糖 尿 病 和 中 风 等 。 例 如 ， 血 管 管 


0 引言 本 
和 径 、 血 管 分 支 角 或 血管 弯曲 度 的 改变 是 高 血压 引起 的 必然 结果 
视网膜 血管 结构 的 变化 是 许多 疾病 出 现 的 临床 表现 ， 这 些 。 .24， 出 现 新 增 血管 是 糖尿 病 视网膜 病变 的 象征 中 。 眼 底 视网膜 


收 稿 日 期 : 2018-01-31; 修 回 日 期 : 2018-03-22 ”基金 项 目 : 国家 自然 科学 基金 资助 项 目 (11501540，11471307); 中 国 科学 院 西部 之 光 基 金 资助 项 目 ; 
重庆 市 教委 科技 项 目 (KJ1501120，KJ1401118) 

作者 简介 : 汪 维 华 (1969-)， 男 ， 重 庆 荣 昌 人 ， 教 授 ， 博 士 研究 生 ， 主 要 研究 方向 为 图 像 处 理 《y2002ww@163.com); 张 景 中 〈1936-)， 河 南 汝 南 人 ， 院 
士 ,博导 ,主要 研究 方向 为 机 器 证 明 、 教 育 数学 、 距 离 几 何 及 动力 系统 ; 吴 文 济 (1976-)， 男 ， 四 川 成 都 人 ,研究员 ,博士 ,主要 研究 方向 为 符号 数值 混合 计算 
.DAE/PDAE. 


201804.02362v1 


chinaXiv 


ChinaxXi V 合 人 期 于 


录用 稿 王 维 等 : 改进 的 形态 学 与 Otsu 相 结 合 的 视网膜 血管 
图 像 在 眼科 疾病 诊断 和 研究 中 得 到 了 广泛 的 应 用 ， 分 割 和 测量 。 然后 提出 一 种 修正 技术 降低 病态 对 血管 的 影响 ， 最 后 使 用 Otsu 
眼底 视网膜 血管 的 结构 状态 在 眼科 疾病 诊断 和 研究 中 提供 了 重 。 方法 对 血管 进行 识别 与 分 割 。 
要 信息 。 因 此 ， 视 网 膜 图 像 血 管 分 割 是 眼科 疾病 诊断 和 分 析 病 > 
情 严 重 程度 的 重要 依据 ， 但 目前 基本 上 还 是 采用 人 工 方式 对 拍 视网膜 血管 分 唱 
摄 的 视网膜 图 像 进行 分 析 ， 这 浪费 了 医生 大 量 的 时 间 ， 分 析 结 。 1.1 图 像 预 处 理 
果 同 医生 的 经 验 相 关 ， 主 观 性 强 。 因 此 ， 视 网 膜 图 像 血管 的 自 一 般 情 况 下 ， 有 眼底 相机 采集 到 的 图 像 是 RGB 格式 的 彩色 
动 分 割 ， 尤 其 是 高 精度 的 自动 分 割 算法 需求 迫切 ， 这 也 是 我 们 。 图 像 ， 这 有 利于 人 工分 割 ， 然 而 对 于 机 器 而 言 ， 大 多 情况 是 处 
的 研究 目的 。 理 单 色 通道 ， 研 究 表明 血管 与 背景 对 比 度 最 强 的 是 视网膜 图 像 
经 过 不 断 探索 ， 人 们 已 经 提出 了 许多 视网膜 血管 自动 分 割 。 的 绿色 通道 即 G 分 量 申 ， 我 们 也 选择 对 G 分 量 进行 研究 处 理 。 
方法 ， 这 些 算法 大 体 上 可 以 分 为 两 大 类 : 有 监督 学 习 的 视网膜 于 图 像 采集 过 程 中 存在 噪声 ， 下 一 步 处 理 之 前 需要 进行 滤波 
分 割 方法 和 无 监督 学 习 的 视网膜 分 割 方法 。 有 监督 学 习 算法 包 ”处 理 ， 主 要 目的 是 消除 孤立 点 噪声 ， 这 里 采用 Gauss 滤波 器 滤 
括 两 大 步骤 ， 训练 阶段 和 识别 分 割 阶段 。 训 练 阶段 利用 人 工分 ，” 波 。 最 后 提取 目标 处 理 兴 趣 区 域 ， 只 对 兴趣 区 域内 部 进行 分 害 
割 好 的 标准 视网膜 血管 图 像 按照 一 定 的 学 习 方 法 训练 分 类 器 。 处理 。 
在 识别 分 割 阶段 使 用 训练 好 的 分 类 器 对 具有 一 定 特征 的 未 知 视 。 | 。 图 像 增强 与 修正 
网 膜 图 像 进行 血管 分 制 。 有 监督 学 习 的 血管 分 割 通常 能 够 获得 Se 甫 是 攻 大 全 纪 
入 下 站 生生 天 度 ， 介 是 该 基 生生 要 证 全 下 在 的 计生， 而 县。 背景 之 癌 的 对 比 度 低 ,这 不 利用 视网膜 图 像 血管 的 分 割 . 因 此 ， 
人 An 和。 在 进行 分 割 之 前 ， 需 要 对 图 像 进行 增强 处 理 。 这 里 采用 改进 的 
鲁 棱 的 混合 特征 集 的 有 监督 学 习 视 网 膜 分 割 算法 ， 算 法 综合 利 。 天 起 志高 低 由 变换 增强 图 像 。 
用 其 他 不 同 算法 的 特征 集合 构成 一 个 混合 特征 向 量 描述 血管 特 
征 ， 由 17 个 特征 构成 17-D 向 量 , 这些 特 征 包括 :13 个 不 同 配 态 学 中 高 帽 变换 (top-hat) 用 开 运 算 对 样本 进行 处 理 , 低 
的 Gabor 滤波 器 、 对 比 度 增强 、 形 态 学 高 帽 变 换 增 强 、 血 管 帆 变 换 (bottom-hat) 则 用 闭 运 算 对 样本 进行 处 理 。 开 运算 能 够 
测量 和 COSFIRE 渡 波 器 。 算 法 使 用 一 个 快速 而 简单 的 随机 森 ”消除 细小 物体 ,平滑 较 大 物体 边界 和 分 离 通过 纤细 连接 的 物体 。 
林 分 类 器 来 训练 和 分 割 血管 。Soares 加 等 提出 了 一 种 基于 Gabor 。 闭 运算 的 作用 是 平滑 边界 、 连 接近 邻 物体 和 填充 物体 等 。 开 运 
小 波 和 贝 叶 斯 分 类 器 的 有 监督 学 习 的 视网膜 血管 分 割 算法 ， 算 算 与 闭 运算 是 形态 学 运算 中 的 基本 操作 局， 分 别 定义 为 
法 使 用 多 尺度 Gabor 小 波 变换 构成 一 个 2-D 特征 向 量 ， 使 用 由 
叶 斯 分 类 器 将 一 个 像素 点 划分 为 血管 和 非 血 管 两 类 。 lm = Lo (1 
无 监督 学 习 的 算法 不 需要 训练 ， 也 不 需要 提供 训练 集 ， 相 
和 上 此 < 从 小 人 Te =7e9。 (O) 
对 而 言 ， 无 监督 学 习 视 网 膜 血 管 分 割 算法 的 执行 速度 快 。 
Hooverl 等 利用 血管 局 部 区 域 信 息 ， 提 出 了 一 种 阔 值 控 针 血管 ”其 中 :表示 图 像 ，$。 表 示 形 态 学 中 的 结构 元 素 。 高 帽 变换 就 是 
检查 方法 ， 在 预 处 理 阶段 ， 算 法 使 用 Gauss 滤波 器 滤波 。 陶 治 。 采用 开 运 算 检查 灰 度 图 像 中 的 灰 度 峰 值 ， 反 之 ， 低 帼 变换 则 是 
江 中 等 提出 了 一 种 基于 相位 一 致 性 的 无 监督 视网膜 图 像 血 管 自 ”采用 闭 运算 检 查 灰 度 图 像 中 的 灰 度 谷 值 趾 ， 高 帽 变换 和 低 帽 变 
动 分 割 算法 ， 该 方法 使 用 改进 的 相位 一 致 性 分 制 血管 ， 再 使 用 换 的 定义 分 别 为 
高 帽 变换 进一步 增强 血管 图 像 ， 最 后 使 用 面积 闵 值 消除 块 状 干 
扰 。 王 晓 红外 等 提出 了 一 种 基于 多 尺度 2-D Gabor 变换 的 视 网 Top pa =1 (To,)=7 -1 G) 
膜 血管 分 割 算法 ， 算 法 使 用 直方 图 均衡 化 进行 预 处 理 ， 然 后 使 
用 多 尺度 Gabor 变换 增强 血管 图 像 ， 最 后 分 别 使 用 形态 学 重 构 Tjonomn pt = — (le) = (4) 
法 和 区 域 增长 法 分 割 粗 细 血 管 ， 将 分 割 的 粗细 血管 融合 获 和 
0 ee 这 里 采用 高 低 帽 变换 组 合 对 图 像 进行 增强 处 理 ， 增 强 处 理 
分 扣 方 法 ， 但 是 这 些 算法 还 是 有 此 地 方 信 每 改进，a) 分 割 的 镶 。 过 各 分 为 两 步 ,第 一 步 分 别 获得 眼底 图 像 的 高 由 变换 和 低 帆 变 
天 性 ， 右 的 算法 分 害 迁 度 科 高 ， 但 是 年 术 性 闫 ， 无 法 笑 应 对 比 。。 痪 图 像 ， 第 二 步 用 高 由 变换 与 低 相 变 换 之 差 以 达到 增强 图 像 的 
四 ee 的 。 图 像 增强 过 程 可 以 表示 为 ， 
度 低 、 病 态 血管 的 分 割 ，b) 分 割 的 精度 低 ， 有 的 算法 分 割 精度 
低 ， 无 法 满足 研究 和 临床 需求 ， 尤 其 是 视网膜 分 支点 和 交点 处 _J _J | 
的 管 分 割 不 饮 准 确 _ enhance top_hat bottom_hat 
针对 以 上 视网膜 血管 分 割 算法 存在 的 问题 ， 提 出 了 一 种 基 了 即 为 增强 后 的 视网膜 图 像 
于 改进 的 基于 高 低 帽 变换 和 Otsu 的 视网膜 血管 自动 分 割 新 算 
法 。 算 法 首先 综合 使 用 高 低 帽 运算 对 眼底 图 像 进行 增强 处 理 ， 


形态 学 变换 中 结构 元 素 的 选择 很 重要 ， 由 于 血管 是 管状 结 


录用 稿 汪 维 华 ， 等 : 改进 的 形态 学 与 Otsu 相 结 合 的 视网膜 血管 分 害 


构 ， 也 就 是 一 种 线性 结构 ， 因 此 ， 运 算 中 的 结构 元 素 5S。 采 用 线 (a) 原始 图 像 (b) 预 处 理 后 的 眼底 图 像 
性 结构 。 该 结构 元 素 有 两 个 参数 : 长 度 (len) 和 角度 〈deg )， 
不 同 参数 的 线性 结构 元 素 模 板 如 图 1 所 示 。 根 据 粗 血管 直径 ， 
本 文中 结构 元 素 的 长 度 取 10 个 像素 点 , 即 len=10。 从 视网膜 图 
像 中 ， 我 们 不 难 发 现 ， 血 管 构成 一 个 网 状 结构 ， 血 管 段 的 方向 
是 不 确定 的 ， 因 此 血管 的 角度 在 0 到 360" 之 间 ， 根 据 对 称 性 ， 


度 选择 为 0 到 170" 之 间 每 隔 10° 选 取 一 个 角度 。 这 样 我 们 可 | 
以 得 到 18 个 模板 进行 高 低 由 变换 结果 , 假设 第 个 模板 的 高 低 II 
由 变换 在 像素 点 Gx) 处 的 灰 度 值 为 ，Joriance ; Ge)， 则 最 终 的 曾 
强 图 像 每 个 像素 点 的 灰 度 值 取 该 18 个 变换 结果 中 对 应 位 置 处 图 2 图 像 的 增强 和 修正 
的 最 大 值 ， 即 1.3 ”Otsu 阅 值 分 割 
本 文 使 用 Otsu 算法 分 割 血管 图 像 。Otsu 算法 [0 也 叫做 最 
Vogt = MaX Lormance it )) (6) 大 类 间 方 差 阐 值 分 割 算 法 ， 是 一 种 动态 国 值 方法 ， 该 算法 最 早 
0 00 0 0 0 0 日 本 学 者 Otsu 于 1979 年 提出 来 的 。 其 基本 原理 是 ， 按照 灰 
000 0 000 1 度 图 像 的 灰 度 值 分 布 ， 将 图 像 划分 为 目标 和 背景 两 类 ， 根 据 两 
000 0 0000 ey 类 之 间 的 方差 进行 分 割 ， 两 类 之 间 的 类 间 方 差 越 大 说 明 这 两 类 
1 1 1 Q@ 111 00 Dm 00 之 间 的 区 分 度 越 大 ， 当 目标 或 背景 被 错 分 时 ， 目 标 和 背景 之 间 
O02 oo 0 1 0 00 的 差别 变 小 , 因此 目标 与 背景 错 分 的 概率 最 小 时 类 间 方 差 最 大 。 
1 ， ， 10 0 0 假设 使 用 闵 值 法 将 一 幅 尺 度 大 小 为 mxn 的 图 像 分 割 为 背景 5 
ay en ld wo a 和 目标 S1 两 类 ， 其 灰 度 级 为 [0L-I] ， 阔 值 为 t+， 则 有 
S60={0,1,2,...，S1={btt1,t+2,.….,L-1} 。 灰 度 级 为 i 的 像素 点 数 
图 1 高 低 由 变换 中 结构 元 素 5 量 表示 为 ns， 则 灰 度 级 i 在 图 像 中 出 现 的 概率 为 
一 个 增强 效果 例子 如 图 2 所 示 ， 图 中 ，(a) 是 一 幅 原始 眼 本 
底 图 像 ，(b) 是 经 过 预 处 理 后 的 视网膜 灰 度 图 像 ，(c) 是 使 用 pi=n,/(mxnm,p; 20, P=1  (®) 
高 低 帆 变换 组 合 增强 后 的 效果 。 从 图 2(c) 可 以 看 出 增强 图 像 。” 则 类 so 和 Si 出 现 的 概率 分 别 为 
中 含有 许多 由 于 病态 引起 的 光照 反射 产生 的 噪声 也 被 同时 增强 
了 。 然 而 ， 仔 细 观 察 原 始 图 像 与 增强 图 像 ， 可 以 发 现在 原始 图 商人 研 > 六 (9) 
像 和 预 处 理 图 像 中 强 光 照 反射 处 背景 与 血管 颜色 对 比 度 非 常 大 ， 
是 不 同 的 颜色 ， 但 是 在 增强 图 像 中 反射 噪声 与 血管 的 对 比 度 却 PD)=D, p;=1- polt) (10) 
非常 小 ， 部 分 强 反 射 光 噪声 也 被 增强 ， 为 此 ， 考 虑 使 用 预 处 理 司 
类 So 和 8) 的 平均 灰 度 级 分 别 为 


上 


图 像 对 增强 图 像 进行 修正 ， 降 低 强 反射 背景 与 血管 之 间 的 对 
度 ， 削 弱 被 误 增强 的 反射 噪声 ， 增 强 血 管 与 背景 的 对 比 度 。 


Re 
WW 


1 /， PP: 
正 的 过 程 如 下 ; Wl) = 2 0 (11) 
0 
1 modify 二 1 preprocess * 1 hat (7) 
ZL-L ,. bp, 
u(t) = > 一 一 ) 12 
其 中 ，Jrevrosece 为 预 处 理 后 的 图 像 ，“.*” 为 图 像 中 的 相应 位 置 之 ! Dm pi(?) 
间 的 位 相 乘 运算 ， 即 数组 对 应 位 置 的 元 素 之 间 的 乘积 。 经 过 修 
正 后 ， 图 2 (ce) 变换 为 图 2 (d)， 由 此 看 出 血管 与 背景 的 对 比 ”“ “人 
度 增加 ， 进 一 步 消 除了 病态 引起 的 高 反射 噪声 干扰 。 a i 


这 样 ， 这 两 类 之 间 的 方差 为 


as 人 = po(ul) uD) + piu) — u(t) 4) 


则 血管 分 割 的 阐 值 为 So 和 Sj 两 类 之 间 方 差 达 到 最 大 值 处 ， 即 
最 优 阔 值 了 为 


T=argmaxfce(D)} (15) 


下 面 分 析 Otsu 算法 的 复杂 度 上 表达 式 可 知 , 在 经 典 的 
Otsu 算法 中 t 的 取 值 范围 为 [1,Z], 假 设 图 像 中 类 So 和 Sj 的 像素 
点 的 个 数 分 别 为 No 和 Nj， 总 像素 点 个 数 为 n 则 对 于 每 一 个 1 
Po(t) 和 uot) 分 别 要 计算 No 次 加 减 运算 和 1 次 除法 运算 , P(t) 和 
u(t) 分 别 要 计算 Ni 次 加 减 运算 和 1 次 除法 运算 , 而 计算 类 间 差 
需要 2 次 乘法 运算 、3 次 加 减 运算 、2 次 平方 运算 。 则 灰 度 级 为 
L 的 图 像 的 运算 为 : Ln+3L 次 加 减 运算 、47 次 乘除 运算 、2L 次 
平方 运算 , 即 算 法 的 操作 次 数 为 T(n)=Ln+9L, 复杂 度 为 O(Ln)。 
为 了 减少 运算 复杂 度 ， 采 用 文献 (中 的 改进 算法 ， 将 t 的 遍历 


范围 缩小 为 [u- a,u+ a]，a=50，a 的 详细 分 析 过 程 参见 文献 021。 


因此 ， 算 法 的 操作 次 数 为 T(n)=50n+450， 复 杂 度 为 O(n)。 
2 ”实验 结果 与 分 析 


2.1 实验 数据 

算法 采用 两 个 公共 视网膜 图 像 数 据 库 DIRVEI9 
进行 测试 DIRVE 图 像 数 据 库 是 由 荷兰 糖尿 病 视网膜 病 
变 组 构建 ， 该 数据 库 一 共有 40 张 图 像 ， 图 像 平分 为 两 个 子 集 : 
训练 子 集 和 测试 子 集 , 每 个 子 集 包 含 20 张 图 片 。 测 试 子 集 包含 
两 个 标准 分 割 数据 库 ， 分 别 由 两 个 专家 手工 分 割 而 成 。 训 练 子 
集 包含 一 个 标准 分 割 数据 库 。 新 算法 是 无 监督 方法 ， 不 需要 训 
练 集 ， 以 测试 子 集 及 其 第 一 个 专家 标准 分 割 库 为 金 标准 进行 血 
管 分 割 和 性 能 测试 与 分 析 。 


与 STARE(3] 


STARE 视网膜 图 像 数 据 库 是 由 Hoover 等 人 采集 整理 而 成 ， 
数据 库 总 共 包含 20 幅 图 像 ， 其 中 10 幅 病理 视网膜 图 像 ，10 幅 
健康 视网膜 图 像 。 数 据 库 包含 两 个 人 工 标注 的 标准 分 制 ， 这 里 
也 是 以 第 一 位 专家 的 分 割 为 金 标准 进行 性 能 测试 与 分 析 。 


2.2 评价 指标 

为 了 对 算法 进行 性 能 测试 ， 就 需要 通过 性 能 指标 客观 地 
行 定量 分 析 。 这 里 采用 特异 性 (specificity， 简 写 为 Spe)、 敏 感 
度 (sensitivity， 简 写 为 Sen) 和 精度 (accuracy， 简 写 为 Acc) 


三 个 指标 进行 分 析 。 这 三 个 指标 分 别 定义 为 "4 站 
Spe=7TN /(FP+TN) (16) 
Sen=7TP/(FN +7P) (17) 
Acc=TN/(IP+FN+rFP+TN) (18) 


其 中 TP (true positive) 为 真 阳 性 ， 表 示 被 正确 分 割 的 血管 
像素 点 个 数 ; TN (true negative) 为 真 阴 性 ， 表 示 被 正确 分 割 的 
非 血管 即 缘 景 像素 点 个 数 ，FP (false positive ) 为 假 阳 性 ， 即 血 
管 被 错误 地 分 割 成 非 血管 的 像素 点 个 数 ; FN (false negative) 
为 假 阴性 ， 即 非 血 管 被 错误 地 划分 为 血管 的 像素 点 个 数 。 
TP+FN+FP+TN 就 是 图 像 中 感 兴趣 区 域 的 总 像素 点 个 数 。 


2.3 ”实验 结果 
为 了 验证 新 算法 的 有 效 性 ,使 | 


MATLAB 14a 平台 对 新 算 
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法 进行 仿真 ， 仿 真 的 机 器 使 用 的 是 主 频 为 2.4 GHz， 
GB 的 PC 机 。 算 法 采用 公开 的 数据 库 DRIVE 和 STARE 进行 
测试 , 图 3 图 4 分 别 展示 了 两 个 数据 库 中 两 幅 视 网 膜 图 像 分 割 
后 的 视觉 效果 ， 图 中 每 一 列 为 一 幅 视 网 膜 图 像 及 其 对 应 的 分 割 
结果 。 3 是 DRIVE 数据 库 的 测试 效果 图 ， 其 中 图 3 第 一 行 
为 原始 视网膜 图 像 ， 图 3 第 二 行为 金 标准 ， 图 3 第 三 行为 新 算 
法 分 割 的 结果 。 图 4 是 STARE 数据 库 的 测试 效果 图 ， 其 中 图 
4 第 一 行为 原始 视网膜 图 像 ， 图 4 第 二 行为 金 标准 ， 图 4 第 三 
行为 新 算法 分 割 的 结果 。 


3 DRIVE 视网膜 图 像 血 管 分 割 效 果 
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2.4 实验 分 析 


图 4 STARE 视网膜 图 像 血 管 分 割 | 部 


常用 的 图 像 分 割 评价 指标 accuracy、specificity 和 sensitivity 对 
岂 分 析 。 在 DRIVE 视网膜 图 


新 算法 进行 定量 


仅 1.5 s。 


节 对 新 算法 与 其 他 算法 进行 比较 与 分 析 。 这 是 
法 与 文献 [18~25] 等 方法 的 性 能 比较 如 表 1 所 示 。 在 STARE 视 
网 膜 图 像 数 据 库 中 ， 新 方法 与 Hoover et al.P9、Jiang et al.081、 


新 算法 的 优势 就 是 对 经 过 ] 
正 ， 修 正 的 了 形态 学 变换 对 视网膜 


日 干扰 ， 图 
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形态 学 变换 的 
入 


像 中 病 


增强 图 像 进行 了 修 
变 对 强 反 射 


le 
Ct 
Ey 


5 为 一 幅 视 网 膜 图 像 修 正 


多。 从 图 5 (c) (Cd 


LH 
DTy 


少 ， 


几 个 


里 选 


昌 库 中 ， 新 方 


割 效果 图 比较 ， 基 中 〈a) 为 原始 视网膜 图 像 ，(\b ) 为 该 图 像 分 
结 


个 
割 的 金 标准 ，(《c) 为 修正 前 分 割 结果 图 ，(d 


前 后 得 到 的 血管 分 


) 为 修正 


修正 后 


) 
\ 体 如 图 中 的 红色 图 内 部 分 所 示 。 


分 割 的 血管 图 中 包 


Lam et al.021、Zhang et al.P0、Martinez-PerezP2U、Bankhead et al.P21、 (c) 修正 前 的 分 割 结果 ”〈d) 修正 后 的 分 割 结 果 
Yin et al.P31]、XueD24、Fraz et al.D51 等 方法 的 性 能 比较 如 表 2 所 图 5 形态 学 变换 修正 的 效果 比较 
示 。 由 表 1 和 表 2 可 以 看 出 ， 新 算法 的 精度 高 于 表 中 其 他 非 监 
督 视网膜 血管 分 割 算法 ， 而 且 新 算法 的 执行 时 间 也 比较 少 ， 仅 为 了 进一步 测试 修正 过 程 对 新 算法 性 能 的 改进 和 影响 ， 我 
们 对 未 添加 修正 过 程 和 添加 了 修正 过 程 进行 了 比较 测试 。 同 时 ， 
四 i ee 我 们 还 进一步 测试 了 多 次 重复 使 用 表达 式 (7) 对 视网膜 血管 图 
表 1 数据 库 DRIVE 中 不 同 视 网 膜 血管 分 割 算法 性 能 比较 像 分 割 效果 的 影响 ， 结 果 如 表 3、4 所 示 。 从 表 3、4 可 以 看 出 
Method Sp Se Time 修正 1 次 的 效果 最 好 ， 修 改 1 次 在 两 个 数据 库 中 都 能 够 得 到 最 
Second observer 0.9470 0.9717 ”0.7796 。” -- 好 的 分 割 精 度 。 同 时 从 表 3、4 中 还 可 以 看 出 ， 修 正 过 程 对 
Jiang et al.[19] 0.9212 ”-- 本 8-36s DRIVE 数据 库 的 影响 较 小 ， 分 割 精度 从 0.9381 变 为 0.9382， 
Lam et al.[20] 0.9579 ”0.7205 。”-- 但 是 对 STARE 数据 库 的 影响 较 大 ， 分 割 精 度 从 0.9403 增加 到 
Zhang et al.[21] 0.9382 ”0.9724 ”0.7120 ”-- 0.9460。 分 析 其 原因 ,STARE 数据 库 中 包含 的 病态 图 像 比 较 多 ， 
Martinez-Perez[22] 0.9344 0.9655 0.7246 -- 而 DRIVE 中 病态 图 像 少 ， 这 也 进一步 表明 新 算法 对 于 人 处理 病 
Bankhead et al.[23] 0.937 0.971 0.703 三 态 图 像 具 有 一 定 优势 ， 也 就 是 新 算法 鲁 棒 性 好 。 
Yin et al.[24] 0.9710 0.6522  -- 表 3 修正 过 程 对 ORIVE 数据 库 的 影响 
Yin et al.[17] 0.9371 。” -- 二 二 
Xu[25] 0.933 0.955 0739  -- Revsion Mune ee 种 Se 
Fraz et al.[26] 0.8112 0.9742 0.7301 -~ 0 :2381 ‘08 038] 
This method 0.9382 0.9926 0.5686 1.6s 1 W9382 00220 9080 
表 2 数据 库 STARE 中 不 同 视网膜 血管 分 割 算法 性 能 比较 ? ed 
3 0.9323 0.9950 0.5069 
Method A wp > Te 4 0.9284 0.9953 0.4739 
PeoRd ODRenYer i 表 4 修正 过 程 对 STARE 数据 库 的 影响 
Hoover et al.gg 0.9275 ”0.81 0.65 Smin 
Jiang et al.[19] 0.9009 ”-- 8-36s Revison numbes Aso Sp Se 
Lam et al.[20] 0.9586 0.6786 -~- 9 20 ln. ed 
Zhang et al.[21] 0.9131 0.9479 0.7716 -- . .3400 "0981S Qe0378 
Martinez-Perez[22] 0.9410 0.9569 0.7506 -- < O00 0 
Bankhead et al.[23] 0.932 0.950 0.758  -- 3 94037 -020 W323 
Yin et al.[24] 0.9666 0.7248 ”-- 4 0 
Xue[25] 0.920 0.944 0.825  -- 3 ”结束 语 
Fraz et al.[26] 0.7294 0.9660 0.7318 -- 
this method 0.9460 0.9815 0.6378 1.6s 本 文 针对 视网膜 图 像 血 管 分 割 中 对 比 度 小 ， 视 网 膜 病 变 在 
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眼底 图 像 采集 过 程 中 反光 干扰 引起 的 光照 反射 过 强 问 题 ， 提 出 
了 一 种 新 的 基于 形态 学 与 Otsu 的 视网膜 血管 分 割 算 法 。 方 法 首 
先 利用 形态 学 中 的 高 低 帽 变换 增强 经 过 预 处 理 后 得 到 的 视网膜 
图 像 , 然 后 提出 了 一 种 修正 增强 图 像 的 方法 解决 强 光 反射 问题 
最 后 采用 Otsu 阔 值 分 割 发 分 割 血管 。 算 法 采用 国际 上 公开 的 
DRIVE 和 STARE 视网膜 图 像 数 据 库 对 新 算法 进行 了 测试 。 实 
验 表 明 新 方法 具有 较 高 的 分 割 精度 ， 能 够 有 效 处 理 视网膜 疾病 
对 血管 分 割 的 影响 , 算法 鲁 棒 性 好 , 算法 的 具有 较 强 的 实用 性 。 
进一步 提高 算法 对 病变 区 域 的 筛 查 以 及 算法 在 视网膜 图 像 分 析 
中 的 应 用 是 我 们 今后 的 研究 方向 。 


Eo 
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